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Abstract. A hydrostatic and stationary white dwarf envelope model has been developed for the study of the 
post-outburst phases of classical novae and their soft X-ray emission. Several white dwarf masses and chemical 
compositions typical for classical novae have been considered. The results show that the luminosity, maximum 
effective temperature and envelope masses depend on the white dwarf mass and on the chemical composition. 
Envelope masses for which equilibrium solutions exist are pretty small (~ 1CP 7 — 1O~ 6 M0), thus leading to a 
short duration of the soft X-ray emitting phase of classical novae, in agreement with most of the observations. 
The models presented provide a useful tool for the determination of the white dwarf properties from observable 
parameters in the X-ray range. 
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1. Introduction 

Classical nova outbursts are caused by the explosive burn- 
ing of hydrogen on the surface of the white dwarf of 
a cataclysmic variable, a close binary system with mass 
transfer from a main sequence star (secondary) to a white 
dwarf (primary). Hydrodynamic models show that, after 
10~ 5 — 10 _4 M Q of H-rich material are transferred from 
the secondary star to the white dwarf, ignition in degen- 
erate conditions takes place in the accreted envelope and 
a thermonuclear runaway is initiated (Starrfield 1989). As 
a consequence of the thermonuclear runaway, the envelope 
expands and causes the brightness of the star to increase 
by ^10 magnitudes, reaching maximum luminosities up 
to ~ 1O 5 L . A fraction of the envelope is ejected at large 
velocities (hundreds or thousands of km/s), while a part of 
it returns to hydrostatic equilibrium with steady nuclear 
burning. Hydrodynamic models also show that a powerful 
outburst is only possible if the accreted matter (with so- 
lar abundances) is enhanced by some mixing mechanism 
with CNO nuclei from the degenerate core. Analyses of 
nebular emission lines from classical novae show indeed 
ejecta enhanced in CNO nuclei compared to cosmic abun- 
dances and, in some cases, enhancement in Ne, Mg and 
other heavier elements. 

The presence of a remaining hydrogen burning enve- 
lope on the post-outburst white dwarf was already pointed 



Send offprint requests to: M. Hernanz 



out after the first hydrodynamic simulations indicated 
that only a fraction of the hydrogen rich accreted mate- 
rial is ejected (Starrfield et al. H974fl . The rest of the enve- 
lope is expected to return to hydrostatic equilibrium with 
steady hydrogen burning, powering a source of constant 
bolometric luminosity. As the ejected envelope expands, 
the photosphere receedes, and inner and hotter layers be- 
come visible, shifting the spectrum from optical to ultravi- 
olet, and finally to soft X-rays. This picture could explain 
the persistence of the high temperature ultraviolet emis- 
sion after the optical decline, detected for the first time 
with OAO-2 in Nova Serpentis 1970 (Gallagher & Code 
ITTJ73T) and in other novae with IUE (1978-1996). 

Without a detailed model for the post-outburst enve- 
lope, the duration of the constant bolometric luminosity 
phase was estimated as the nuclear timescale of the mass 
left on the white dwarf surface (see for instance table 1 in 
Gchrz ct al. 1998). Since the outburst ejects a fraction of 
the accreted envelope, the remaining envelope mass was 
assumed to be of the same order of magnitude, which was 
estimated as the envelope mass needed to achieve the crit- 
ical pressure at the base of the envelope to trigger the 
outburst (MacDonaldHnSI, ~ 10~ 4 - 10" 5 M Q for white 
dwarf masses from 0.9 to 1.3 M Q (in agreement with ac- 
creted envelope masses obtained with hydrodynamic mod- 
els). With the expected constant bolometric luminosity, 
10 4 — 10 5 L Q , the nuclear time-scale during which the nova 
was expected to be a bright soft X-ray source was of tens or 
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hundreds of years. According to these estimations, classi- 
cal novae would be expected to shine as soft X-ray sources 
for decades or centuries. The spectrum of this super-soft 
X-ray emission is expected to be that of a hot white dwarf 
atmosphere (MacDonald & Vennes l*L991[) . 

The first detection of soft X-rays from a classical nova 
in outburst came from GQ Mus 1983, when observed 
with EXOSAT in 1984, 460 days after optical maximum 
(Ogelman et al. IT384). Other novae were also detected 
some hundred days after optical maximum (PW Vul, QU 
Vul, Ogelman et al. 11984(1 . but a major break-through was 
provided with the launch of the Rontgcnsatellit (ROSAT), 
which observed a total of 39 recent classical novae, less 
than 10 years after the outburst. A systematic search for 
the X-ray emission from those post-outburst classical no- 
vae in the ROSAT archival data was performed by Orio et 
al. I|2001[) . Contrary to what was expected, very few no- 
vae were discovered as X-ray emitters shortly after their 
outburst. In fact only 3 novae showed a soft X-ray spec- 
trum: GQ Mus, V1974 Cyg and Nova LMC 1995. GQ Mus 
showed the longest constant bolometric luminosity phase 
so far observed: its soft X-ray emission lasted for 9 years, 
still less than predicted. V1974 Cyg, the best monitored 
soft X-ray nova, turned off only 18 months after the out- 
burst. In the case of Nova LMC 1995, it was detected by 
ROSAT as a soft X-ray source three years after outburst, 
and in a second observation with XMM-Newton in the 
year 2000, five years after outburst. 

X-ray observations give direct insights into the hydro- 
gen burning layer on the white dwarf surface after the 
outburst, but the presence of the soft X-ray source can 
also be determined from ultraviolet observations (Shore 
Cassatella et al. EStMjl . The ultraviolet is the key 
spectral regime for studying the dynamics and opacity of 
the expanding shell, ionized by the central X-ray source, 
and determine abundances, velocities, spatial structure 
and mass of the ejecta. The ultraviolet emission lines are 
indicators of the ionization and density structure of the 
ejecta, and their variations reflect the evolution of the ex- 
panding shell, but also that of the luminosity and effective 
temperature of the central ionizing X-ray source. This in- 
direct indicator has been used to determine the turn-off of 
classical novae from IUE observations (Shore et al.[l996 
Gonzalez- Riestra et al. 119981 Vanlandingham et al. I2001|l . 
In all cases, turn-off times were smaller than expected, 
and also smaller than the 9 years of GQ Mus. 

The discrepancy between theory and observations 
pointed out that an extra mass loss mechanism should 
be present shortly after the nova outburst. Some mech- 
anisms have been suggested, involving the interaction of 
the ejecta with the envelope of the secondary and/or thick 
winds. MacDonald et al. (1985 ) calculated that a fraction 
of the envelope could be lost by dynamical friction of the 
secondary in a common envelope phase following the nova 
outburst. Kato & Hachisu l|1994f) calculated a sequence 
of optically thick wind solutions that simulated the decay 
phase of the visual light curve of novae. 



Nevertheless, the discrepancy between observations 
and theory arose from a comparison mainly based on the 
turn-off times. A model for post-outburst envelopes to 
which X-ray observations could be compared was lack- 
ing. However, a good number of works were devoted to 
the study of hot white dwarf envelopes. Stellar evolution 
models following the post AGB phase simulate a similar 
situation. For population I stars (with Xh=0.7), Paczyhski 
(1970) showed that after exhaustion of He in the core, 
stars of several initial masses (between 3 and 7 M Q ) fol- 
low a common track in the HR diagram with constant lu- 
minosity. Effective temperature increases along the track, 
while hydrogen and helium shell burning reduce the en- 
velope mass. Paczyhski (1970| was also the first to find 
that luminosity is directly proportional to the core mass 
during this phase, that lasts as long as the envelope mass 
is larger than some critical value. When this critical point 
is reached, the shell sources die out and the degenerate 
core starts to cool down. 

A number of studies have been devoted to accret- 
ing white dwarfs and the development of thermal pulses 
or flashes on a hydrogen rich envelope. Among them, 
Iben ijH^}) and Iben & Tutukov ijTi^jl found a core 
mass - luminosity relation for accreting white dwarfs with 
steady hydrogen burning in the accreted envelope (with 
Xh=0.64) similar to Paczyhski's one. Iben l)1982|l sim- 
ulated the evolution of a white dwarf envelope without 
accretion (or very small values of accretion rate) as a se- 
quence of steady-state models, and found that the con- 
stant luminosity phase lasted for decades, in agreement 
with the nuclear timescales expected for the envelope 
masses of post-outburst classical novae predicted by hy- 
drodynamic models. The steady hydrogen burning solu- 
tions were combined with an evolving quasi-static approx- 
imation to study the thermal pulses occurring on accret- 
ing white dwarfs. Different points of view of the problem 
can be found in Fujimoto l|1982[) and Paczyhski (1983), 
who developed simple approximations for the study of 
the hydrogen burning envelopes stability, and the evolu- 
tion of shell flashes for different accretion rates. Taking 
into account the existing numerical models for degener- 
ate stars with hydrogen-burning shells, Iben & Tutukov 
(1996) pointed out that the core mass - luminosity rela- 
tion was not unique, and that it depended on the thermal 
history of the star. They compared the models with data 
from symbiotic stars and cataclysmic variables, exploring 
the effects of mechanisms not present in the models (such 
as wind mass loss and mixing mechanisms) in the evolu- 
tion of accreting white dwarfs. Nevertheless, none of those 
works studied envelopes with abundances typical for clas- 
sical novae. Tuchman & Truran (1998) studied the com- 
position influence on the core mass - luminosity relation 
and its implications for the turn-off times of classical no- 
vae. They applied their relation to GQ Mus and V1974 
Cyg, and concluded that the envelope masses required to 
initiate the thermonuclear runaway were larger than the 
envelope masses below which stable models for steady hy- 
drogen burning exist. This condition led them to suggest 
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that the remnant envelope masses in the post-outburst 
stage could be significantly reduced by some kind of dy- 
namical instability. 

In this paper, a numerical model for white dwarf en- 
velopes with steady hydrogen burning is presented. The 
model provides a composition dependent core mass - lu- 
minosity relation, and is mainly aimed to simulate the 
conditions of the remaining hydrogen-rich material left 
on the white dwarf after the nova outburst. Photospheric 
properties are obtained for comparison with X-ray obser- 
vations, and their evolution is approximated as a sequence 
of steady-states. Section 2 describes the envelope model, 
while the main results are presented in section 3. In section 
4, analytical relations between the main global properties 
of the envelopes are given. The quasi-static evolution is 
presented in section 5, and a summary and discussion of 
the results can be found in section 6. 

2. The white dwarf envelope model 

A grid of white dwarf envelope models with steady hy- 
drogen burning has been computed for five different core 
masses (from 0.9 to 1.3 M Q ) and four different chemical 
compositions, all of them corresponding to realistic en- 
velopes accreted onto white dwarfs in novae. Compositions 
are taken from Jose & Hernanz l|1998[) . who developed hy- 
drodynamic simulations of nova explosions with detailed 
computation of the nucleosynthesis, and obtained results 
for the composition of the ejecta compatible with ob- 
served abundances in many cases. As already found from 
the first hydrodynamic models of nova outbursts, their 
models showed that the material accreted from the sec- 
ondary, with solar abundances, needs to be enhanced with 
CNO nuclei by some mixing mechanism between the en- 
velope and the degenerate core. The compositions consid- 
ered here correspond to CO white dwarf envelopes with 
50% mixing of the solar accreted matter with the core 
material, and to ONe white dwarf envelopes with 25%, 
50% and 75% mixing. For the calculation of the accreted 
envelope composition resulting from the mixing, Jose & 
Hernanz l|1998J) took the composition of the degenerate 
cores from Ritossa et al. (1996) for the ONe white dwarfs, 
and assumed 49.5% of oxygen, 49.5% of carbon and 1% of 
neon-22 in the CO cases. Table ri] shows the details of the 
chemical compositions. 

2.1. Numerical Method and Boundary Conditions 

The numerical code integrates the equations of stellar 
equilibrium from the photosphere to the base of the en- 
velope. The envelope is tested at each step for convec- 
tive instability according to the Schwarzschild criterion. 
When convection appears, the temperature gradient is 
solved using the standard mixing-length theory of convec- 
tion (B6hm-VitenseriJ)58). Radiative opacities have been 
calculated for each composition with the on-line version 
of the OPAL code from the Lawrence Livermore National 
Laboratory (Iglesias & Rogers IT99ri|l . and the energy gen- 



eration rates for hydrogen burning are the analytcal ap- 
proximations from Reeves l|1965fl and Paczyhski (1983). 
The equation of state is that of an ideal gas with radiation 
pressure. It has been checked that for the gas conditions 
in the present models, no significant degeneracy (at the 
base of the envelope) or partial ionization effects (at the 
outermost layers) need to be taken into account, and that 
the energy production due to helium burning is negligible. 

The stellar equilibrium equations are integrated in- 
wards using an embedded (adaptative stepsize controlled) 
Runge-Kutta method with Cash-Karp parameters (Cash 
& Karp 119901 Press et al.Q.992). Integration starts at the 
photosphere, with a fixed value for the total luminosity 
and trial values for the total mass of the envelope and 
the photospheric radius, and proceeds inwards until the 
luminosity reaches a sufficiently small value, which is con- 
sidered to indicate the base of the hydrogen burning layer. 
As will be seen later, for all white dwarf masses and com- 
positions there is a plateau of constant luminosity where 
it experiences little change for a wide range of effective 
temperatures. In this plateau, the total luminosity does 
not determine a unique solution anymore and the photo- 
spheric radius is fixed instead. 

The boundary conditions at the photosphere are ob- 
tained from a grey atmosphere model of the same com- 
position as the envelope, which provides the value of the 
atmosphere pressure at the photosphere for the fixed value 
of the luminosity and the trial values of the envelope mass 
and photospheric radius. The atmosphere is tested for con- 
vective instability during integration, which has never ap- 
peared before reaching the photosphere in this work. At 
the base of the envelope, the envelope mass is required 
to be smaller than 10 _13 M Q , while the radius, Rb aS o, is 
compared to the core radius, R c , by requiring the factor 
Rbaac-Rc to j-, e gma n er than 10 -5 . If these conditions 

are not fullfiled, the integration restarts from the photo- 
sphere with new trial values for the envelope mass and 
photospheric radius. The core radius is taken from the 
Hamada & Salpeter l|1961f) core mass-radius relation for 
zero temperature white dwarfs. We have checked that a 
slight increase in the core radius does not affect our main 
results; this radius increase is only expected if degeneracy 
is lifted in the outer core shells as a consequence of heat- 
ing, which is not probable since the timescale for heating 
propagation towards the core is quite large (Henyey & 
L'Ecuyer. rjnSjl . 

3. Results 

The luminosity versus the effective temperature for all 
models, i.e. their position in the HR diagram, is shown 
in left panel of figure ^ For each composition and white 
dwarf mass, the envelope solutions have a maximum ef- 
fective temperature that divides the series of models into 
two branches: a branch of high, quasi-constant luminosity, 
where the maximum value for the luminosity is reached; 
and a branch of low luminosity and constant photospheric 
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Table 1. Chemical compositions of the white dwarf envelope models, in mass fractions. Z s contains metals (all 
elements except H and He) in solar fraction. SXi indicates extra mass fraction of element i beyond that included in 
Z s . The total abundance of metals is indicated by Z tota i. 
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Fig. 1. Total luminosity (left panel) and envelope mass (right panel) versus effective temperature for our white dwarf 
envelope models: ONe75 (dotted line), ONe50 (short dash), ONe25 (long dash) and CO50 (short dash - dot). Effective 
temperature is given in eV, kT c g(eV) = 8.617 x 10~ 5 T o ff (K) Solid lines indicate envelopes without convective regions. 
Five series of models are plotted for each chemical composition, corresponding to total masses 0.9, 1.0, 1.1, 1.2 and 1.3 
Mq . A line indicating the luminosity-effective temperature relation for constant photospheric radius is over-plotted in 
left panel for comparison. 



radius. The models of the high luminosity branch have in 
all cases a convective region below the photosphere, while 
the energy transfer in most of the envelopes of the low lu- 
minosity branch is radiative in the whole structure. Right 
panel in figure ^ shows the envelope mass as a function of 
effective temperature. For every white dwarf mass there is 
a minimum value for the envelope mass, which occurs in 
the vicinity of the maximum effective temperature. The 
envelope mass along the high luminosity branch is almost 
constant, decreasing only slowly for increasing effective 
temperatures, while the decrease of the envelope mass for 
increasing effective temperatures along the low luminosity 
branch is steeper. A summary of the main properties of 
all model envelopes is found in table [5] 



3.1. Stability 

A stability analysis of the shell source indicates that the 
burning shell is only thermally stable for the solutions 
along the high luminosity branch. The unstability of the 
low luminosity branch can also be understood in a quali- 
tative way if the evolution of the envelope is approximated 
as a succession of steady states with decreasing envelope 
masses. Figure |21 shows the temperature at the base of the 
envelope for CO50 models. The high luminosity branch 
in figure ^ corresponds to envelopes with high base tem- 
peratures in figure |21 while envelopes with low base tem- 
peratures correspond to the low luminosity branch. From 
left panel in figure Q]it is clear that, along any of the two 
branches, the envelope mass decreases towards higher ef- 
fective temperatures. If the evolution of an envelope is ap- 



G. Sala and M. Hernanz: Envelope models for post-outburst novae 



5 




10 -7 I ' 1 1 1 1 

4xl0 7 6xl0 7 8xl0 7 10 s 

T base (K) 

Fig. 2. Envelope mass versus base temperature for CO50 
models. Solid lines indicate envelopes without convective 
regions. 



proximated as a succession of steady states, it will proceed 
in any case towards decreasing envelope masses, i.e., in- 
creasing the effective temperature. Hydrogen burning will 
continue until the minimum envelope mass and maximum 
effective temperature are reached; evolution can not pro- 
ceed further with stationary hydrogen burning, since no 
equilibrium configuration for a lighter envelope exists. 

Along the high luminosity branch, an envelope evolves 
at quasi-constant luminosity, increasing the effective tem- 
perature as a consequence of the decrease of the photo- 
spheric radius, which sinks the photosphere deeper into 
the envelope. Nevertheless, the envelope as a whole is 
slowly cooling down and the energy production decreases 
with time, with both the luminosity and the temperature 
at the base of the envelope slowly decreasing (figures Q 
-left panel- and[5J|. But along the low luminosity branch 
(the low bottom temperature branch in figure 0), for de- 
creasing envelope masses, the temperature at the base and 
the luminosity of the envelope would increase. This would 
require an extra source of energy and indicates that evo- 
lution cannot proceed along this branch. 

The stable, high luminosity branch is limited at its low 
effective temperature end by the Eddington luminosity. 
For the low temperature models along this branch, the ra- 
diative luminosity closely approaches the Eddington limit 
as effective temperature decreases, mainly because there's 
a photospheric opacity increase leading to an Eddington 
luminosity decrease. Therefore, envelopes in hydrostatic 
equilibrium do not exist for effective temperatures lower, 
and thus envelope masses and photospheric radii larger, 
than those shown in figure ^ 



3.2. Interior structure 

Figure [3] shows the interior structure of some of the 
CO50 and ONe50 envelope models of the high luminosity 
branch, all of them with the same effective temperature 
(~ 4.6 x 10 5 K ~ 40 eV). The temperature and luminosity 
profiles show that only the deepest layers have tempera- 
tures high enough to contribute to the energy production. 
Luminosity stays constant in most of the envelope. The 
highest temperatures at the base are needed for the most 
luminous envelopes, and in a similar way the most massive 
envelopes produce the highest pressures at the bottom. 

As mentioned above, envelopes on the high luminos- 
ity branch have a convective region just below the pho- 
tosphere, where opacity reaches its largest values. Since 
convection transfers energy outwards more efficiently than 
radiation, the temperature gradient is flatter in the con- 
vective regions, as well as pressure and density gradients, 
contributing to a larger photospheric radius and luminos- 
ity. It can be observed that for the largest luminosities 
and most massive white dwarfs, the density gradient in 
the convective region is close to zero. For high luminosi- 
ties and gravities, inversions of the density gradient have 
been frequently found in convective layers of white dwarf 
envelopes when using the mixing length theory (Fontaine 
& van HornHnH Henyey et al. H5S5l Latour lHTTUl Mihalas 
1965). However, although the density gradient is very flat 
in some of the present models, such inversions do never ap- 
pear. The convective region is more extended for increas- 
ing white dwarf masses. As a consequence, temperature, 
pressure and density gradients in the inner radiative layers 
are steeper for massive white dwarfs, causing the energy 
production layer to be thinner in these cases, which is 
clearly visible in the luminosity profiles. In a similar way, 
the envelope mass increases only in the deepest layers in 
massive white dwarfs, since most of the envelope has small 
density values. 

The effects of the different chemical compositions on 
the interior structure are shown in figure 01 where profiles 
for five models with different abundances but the same 
core mass and total luminosity are plotted. The strongest 
difference in behaviour corresponds to the hydrogen poor- 
est model, ONe75. In this case, the larger metalicity cause 
the opacities to be larger than in hydrogen richer en- 
velopes, causing the convective region to be much more 
extended. As a consequence, profiles in the radiative re- 
gions are steeper, specially in the deepest layers. 

Different chemical compositions and white dwarf 
masses combinations give envelopes with the same observ- 
able properties (luminosity and effective temperature), 
as can be seen in figure Q Figure El plots the interior 
structure of three ONe envelopes with the same effective 
temperature and luminosity ( ~ 5.2 x 10 5 K ~ 45 eV, 
~ 5.2 x 10 4 L Q ), and thus the same photospheric radius 
(2x 10 9 cm), but different white dwarf masses and composi- 
tions (ONe25 with 1.3M Q , ONe50 with 1.2M Q and ONe75 
with 1.1M Q ). In spite of the differences in abundances, the 
profiles are very similar, resulting in equal photospheric 
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properties. The convection region is more extended for 
the hydrogen richest model, ONe25, which is in general 
cooler (except for the innermost layers). For comparison, 
an ONe50 model of the unstable, low luminosity branch 
and with the same effective temperature is also plotted. 

There are small differences in the general properties of 
the CO50 and the ONe50 envelope models, despite hav- 
ing the same hydrogen mass fraction. From figure Q it is 
clear that, for the same effective temperature, the CO50 
models have slightly higher luminosities and larger enve- 
lope masses. Equivalently, for the same luminosity, CO50 
envelope models have higher effective temperatures and 
smaller photospheric radii. This is evident in figure [21 
where CO50 and ONe50 models are plotted for several 
masses, all of them with the same effective temperature. 
For the same effective temperature, distributions of tem- 
perature, pressure and density (in fraction of total enve- 
lope radius) are steeper for ONe50 models. This is directly 
related to larger opacities for ONe50 than for CO50 com- 
positions in the temperature range of the radiative region 
of the envelopes, with T> 10 6 K, that cause the tempera- 
ture distribution to change more rapidly with radius. Since 
temperature decreases more rapidly in the inner regions 
of ONe50 models than for CO50, the energy production 
layer is thinner and the total luminosity of the envelope 
is smaller. Analogously, the density profile in the inner 
regions is also steeper for ONe50 models, causing the en- 
velope mass to increase only in a narrower region than in 
CO50 models and the total envelope mass to be smaller. 

To check that the differences between CO50 and 
ONe50 models are due to the differences in the opacities, 
some envelope models have been reintegrated with opaci- 
ties artificially altered by a certain factor. We have found 
that, as in the case of the CO50-ONe50 comparison, en- 
velopes with higher opacities have smaller luminosities and 
envelope masses. A global increase (decrease) of the opac- 
ity in the envelope by a 10% causes a decrease (increase) 
in the envelope mass of 4% and a decrease (increase) in 
the luminosity of 10%. 

4. Analytical relations between main envelope 
properties 

4.1. Core mass - luminosity relation 

From figure Q it is clear that luminosity is almost con- 
stant along the stable, high-luminosity branch, and that 
its value increases with the core mass. In the left panel 
of figure El plateau luminosities are plotted against core 
mass for all models. The resulting core mass - luminosity 
relation for our ONe models can be approximated by 

£oNc CaU (£o) - 5-95 x 10 4 (J^- - 0.536X H - 0.14) (1) 
and for the CO50 models by 

L p £T(Lo) - 5.95 x 10 4 (j^ - 0. 3 ) (2) 



The plateau luminosity increases with the white dwarf 
mass, and for increasing mixing factors of the accreted 
solar matter with the core material (i.e., for decreasing 
hydrogen abundances). Envelopes with a lower hydrogen 
mass fraction need higher temperatures to reach an equi- 
librium configuration, and therefore have larger luminosi- 
ties. The plateau luminosity versus the hydrogen mass 
fraction is shown in the right panel in figure|Blfor all mod- 
els. Relation is over-plotted for each core mass. While 
ONe models follow (£Q) (except for the most massive white 
dwarf), CO50 models follow a parallel relation but are 
slightly more luminous than ONe50 (between a 3% and a 
7% more luminous), which can also be seen in figure 

The linear relation between the luminosity of a shell 
burning star on a degenerate core and its core mass has 
been known since long. An analytical justification for this 
linear dependence was obtained by Kippenhahn (1981) 
(see also Kippenhahn & Weigert Il990|) . He showed that 

the luminosity was directly proportional to the core mass 

p 

and independent of its radius for p <C 1 , where p = ■ 
With this analytical work, Kippenhahn 1|1981[) confirmed 
the results found with numerical evolutionary calculations 
by Paczyhski l)197U|l . who showed that the luminosity of 
post-AGB stars with a double burning shell (H, He) on 
top of a mostly degenerate CO core could be expressed as 
L{L Q ) = a(M c /M Q - b), with a = 59,250 and b = 0.522. 
Later works found similar results, with small changes in 
the parameters. Among them, Iben l|1982[) . found that in 
the case of steady state, a = 59,500 and b = 0.5. For the 
hydrogen mass fraction used by Iben 1|1982|) (Xh=0.64), 
the 0.536 Xh+0.14 factor in equation Q), 0.48, is very 
similar to the value found in that work, 0.5. 

The increase of the luminosity with decreasing hydro- 
gen mass fraction was also found in analytical studies by 
Tuchman et al. (1983) and with numerical calculations of 
steady hydrogen burning on white dwarfs by Tuchman & 
Truran (1998). They studied the dependence of the lumi- 
nosity with the hydrogen mass fraction for a fixed metal- 
icity (Z=0.25), and for a constant He/H ratio (assumed 
to be solar, « 0.1), and summarized their dependences 
by the relation L(L Q )=5.2xl0 4 [M c /M Q -0.205-0.5(X-Z)]. 
According to the authors, this relation holds also for other 
metalicities. More recently, Marigo H2000f) used an homol- 
ogy formalism to justify the linear core mass - luminosity 
relation and study the dependence of the luminosity with 
the chemical composition, for giant stars with core masses 
between 0.5 and 0.7 M . 

4.2. Other relations 

The maximum effective temperature also increases with 
decreasing hydrogen mass fractions and with increasing 
core masses (see figured), but in this case, no clear relation 
has been found. As shown in the figure, the maximum 
effective temperature does not follow any unique linear 
relation or power law with the core mass. Its behaviour can 
be only approximated by a power law, fcT e ff(eV) ex M^ 10 
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(for core masses smaller than l.lM©) and kT D g(eV) oc 
M 5 J 9 (for more massive white dwarfs). 

As in the case of the luminosity, the envelope mass in 
the luminosity plateau depends on the core mass and on 
the hydrogen mass fraction (see figure^ right panel). Iben 
( 1982} found that the relationship between the envelope 
mass and the white dwarf mass in the luminosity plateau 
was logM££ nS2 (M®) ~ -4.38 - 4.96(M C /M - 1). Since 
the envelope mass decreases along the luminosity plateau 
for increasing effective temperatures (see figure ^) a clear 
definition of the envelope mass can not be established. In 
figure |H1 (left panel), the maximum envelope mass for our 
stable models is plotted as a function of the core mass. It is 
clear that for our models, the logM^ x — M c relation is not 
linear. A general relation between the maximum envelope 
mass, the core mass and the hydrogen mass fraction can 
be found for our envelopes: 

logM-f(M )~O.42X H - (|^-0.13) -5.26 (3) 

which is indicated in figure |S] 

The minimum envelope mass does not have a strong 
dependence on the hydrogen mass fraction (see right panel 
of figure |8J), and can be approximated by the relation 

log M™»{M Q ) ~ 0.13X H - (j^ - 0.13^ - 5.26 (4) 



5. Quasi-static temporal evolution 

An estimation of the evolution of a steady hydrogen burn- 
ing white dwarf envelope in the plateau of quasi-constant 
luminosity can be approximated by a series of station- 
ary envelope models with decreasing envelope masses. The 
layers of the envelopes whose hydrogen has been processed 
are assimilated to the core and do not take part in the in- 
tegration of the next envelope model (the increment in 
the core mass is negligible). The time At needed for the 
envelope to change its mass AM env only due to hydrogen 
burning can be estimated as 



At 



AAf cnv XH 



(5) 



where Xh is the hydrogen mass fraction, L is the mean lu- 
minosity in this period (which varies little in the constant 
bolometric luminosity phase) and e is the energy released 
per gram of hydrogen processed (e = 5.98 x 10 18 erg g _1 ). 
In figures IIUIIII time intervals needed for the envelope 
to evolve between two adjacent marked points are indi- 
cated in days. Evolution proceeds slower for smaller white 
dwarfs, with lower luminosities. Time intervals are also 
longer for envelopes with increasing hydrogen abundance, 
which have smaller luminosities. 

The times found in this work are much shorter than 
the evolutionary times usually found in the literature for 
stationary hydrogen burning on white dwarfs, including 



for instance the times found by Iben (1982). The differ- 
ence can be explained by the smaller hydrogen fraction 
in the envelope: the equilibrium configuration is hotter 
and brighter for hydrogen poorer envelopes, and envelopes 
masses are smaller. To keep this higher luminosity, hy- 
drogen is consumed faster and, in addition, for the same 
amount of hydrogen consumed, a larger fraction of the en- 
velope is reduced. As envelope masses are smaller and are 
reduced at a faster rate, the total timescale to consume 
the envelope is shorter. 

These timescales can be compared with the nuclear- 
burning lifetime for classical novae hydrostatic remnants 
found in the literature. In Starrfield (JTJ2H2I, this lifetime 
was estimated as 



= 400 



10" 4 M Q 



L 



2 x 1O 4 L 



(6) 



where Mh is the hydrogen mass and L the total lu- 
minosity. A similar expression, only a factor | larger, 
was given by MacDonald et al. (1985). In both cases, 
the nuclear lifetime for an envelope of 10~ 5 — 10~ 4 M Q 
is of decades or even centuries. According to hydrody- 
namical simulations, 10~ 5 — 10~ 4 Mq is the mass required 
for the thermonuclear runaway to be triggered (Jose & 
Hernanz \MB Starrfield IT5M1 Prialnik & Kovetz 
and the fraction of the accreted mass that is expected to 
be ejected is not large enough as to reduce the mass of the 
remnant by two orders of magnitude. But in the present 
work, stable envelopes have been found to have masses of 
~ 10~ 6 Mq, which requires an extra mass loss mechanism 
after the nova outburst for the envelope to reach the equi- 
librium configuration. Once the envelope has reached the 
steady hydrogen burning phase, the nuclear lifetime can 
be rewritten with our typical values (M cnv = 1O _6 M0, 
L=4xl0 4 L m ) as 



= 2Xv 



M c . 



L 



4 x 10 4 L C . 



yr 



(7) 



For X=0.35 (as in ONe50 and CO50 models), r nuc =0.7 
yr~ 255 days, which agrees with the times in figure ITUI 

An indication of the typical evolutionary times for our 
models is given in the last column in tabled which lists 
the time needed for the envelope to change its effective 
temperature the last ~ 10 eV before reaching the maxi- 
mum effective temperature. 

6. Summary 

A numerical model has been developed to simulate the 
physical conditions in the steady hydrogen burning enve- 
lope of a white dwarf after a nova outburst. A grid of en- 
velope models has been computed for five different white 
dwarf masses, from 0.9 to 1.3 M Q , and four chemical com- 
positions, corresponding to realistic CO and ONe white 
dwarf envelopes, with several degrees of mixing between 
the core material and the solar-like accreted matter. The 
results show that there exists a maximum luminosity and 
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Table 2. Main properties of the envelope models. 



Model 


M c 


T plclt C3.L1 


1 j | ■ y\\ ax 

ki eff 


rpmax 


TV r Til a V 

M 

lvl cnv 


Tyrmm 
ivJ -cnv 


A t 10,eV 




(Mr,) 




feV) 


C10 5 K ) 


ClO~ 6 Mn) 


flO~ 6 Mn1 


( davsl 




0.9 


3.9 


57 


6.6 


2.3 


1.9 


47 


ONe75 


1.0 


4.4 


64 


7.4 


1.4 


1.2 


24 


(X=0.18) 


1.1 


5.2 


73 


8.5 


0.84 


0.71 


12 




1.2 


5.7 


86 


9.9 


0.40 


0.35 


4.6 




1.3 


6.2 


103 


11.9 


0.16 


0.15 


1.1 




0.9 


3.5 


53 


6.1 


2.7 


1.5 


160 


ONe50 


1.0 


4.1 


60 


7.0 


1.7 


1.2 


78 


(X=0.35) 


1.1 


4.7 


69 


8.0 


0.99 


0.72 


37 




1.2 


5.2 


81 


9.4 


0.45 


0.35 


14 




1.3 


5.6 


98 


11.4 


0.19 


0.15 


4.9 




0.9 


2.9 


49 


5.7 


3.0 


2.2 


430 


ONe25 


1.0 


3.5 


56 


6.5 


2.1 


1.4 


210 


(X=0.53) 


1.1 


4.0 


64 


7.5 


1.20 


0.81 


98 




1.2 


4.6 


77 


8.9 


0.60 


0.39 


36 




1.3 


5.2 


93 


10.8 


0.22 


0.16 


12 




0.9 


3.6 


54 


6.2 


2.9 


1.9 


230 


CO50 


1.0 


4.2 


61 


7.1 


1.8 


1.2 


90 


(X=0.35) 


1.1 


4.7 


70 


8.1 


1.0 


0.72 


48 




1.2 


5.4 


82 


9.5 


0.49 


0.35 


19 




1.3 


6.0 


99 


14.7 


0.19 


0.15 


6 



Time needed for the envelope to evolve from kT e ff 



lOeV to kTX 



a maximum effective temperature for every white dwarf 
mass and composition. The maximum effective tempera- 
ture divides the series of models into two branches: a sta- 
ble branch of high, quasi-constant luminosity, and a low- 
luminosity unstable branch, with constant photospheric 
radius. For every white dwarf mass and composition there 
is also a minimum envelope mass, which occurs in the 
vicinity of the maximum effective temperature. The mass 
of the stable branch envelopes is almost constant, decreas- 
ing only slowly for increasing effective temperatures. The 
plateau luminosity, maximum effective temperature, and 
maximum and minimum envelope masses for each model 
are summarized in tabled as well as the typical evolution- 
ary timescales. Evolutionary timescales are longer for less 
massive white dwarfs, since luminosity is lower in these 
cases. The evolution also proceeds slower for envelopes 
with increasing hydrogen abundances, for the same rea- 
son. 

Each of the overall properties of the envelopes is re- 
lated to the total white dwarf mass and its chemical com- 
position, and we have found analytical approximations 
for our ONe envelope models (equations 1114(1 . Our results 
show that masses of stable envelopes with steady hydro- 
gen burning (~ 10~ 6 M Q ) are in general smaller by at least 
one order of magnitude than envelope masses needed to 
trigger the outburst (~ 1O _5 M0), according to hydrody- 
namical models. Since those models predict that only a 
fraction of the accreted envelope is ejected, it is clear that 
some mass loss mechanism must act after the outburst, 



before the steady hydrogen burning can be established in 
the envelope. An interesting possibility for the mass loss 
is, as pointed out by Tuchman & Truran l(1998J) . the dy- 
namical instabilities of the envelope itself when it is left 
after the outburst with a mass much larger than that of 
stable configurations. 

With the masses found for stable envelopes, the nu- 
clear time-scales are of some weeks or months, similar to 
observed for post-outburst classical novae. For accreted 
envelopes highly mixed with the degenerate core, the short 
duration of the steady hydrogen burning phase could ex- 
plain the lack of soft X-ray emission in most novae ob- 
served, even if not all the accreted envelope is ejected by 
the outburst. For those cases where the soft X-ray emis- 
sion has been detected, the photospheric properties ob- 
tained with the envelope models and their temporal evo- 
lution can be directly compared with the observational 
results. The comparison of observations with the present 
models has been performed in Sala & Hernanz (2005) for 
the ROSAT observations of V1974 Cyg. From this com- 
parison, the white dwarf mass, and envelope mass and 
composition have been constrained. 

The present envelope models represent a powerful tool 
for the study of classical novae, specially if their post- 
outburst phases are monitored in X-rays. With short and 
frequent observations, with a good sensitivity in the soft 
X-ray range and a moderate spectral resolution, the com- 
parison with our envelope models makes possible the de- 
termination of the main parameters of the white dwarf. 
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Our envelope models can also be applied to the study of 
other systems with stationary hydrogen burning on top of 
a white dwarf, (such as supersoft X-ray sources or symbi- 
otic stars). 

Acknowledgements. This research has been partially funded by 
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Fig. 3. Interior structure of ONe50 (dashed line) and CO50 (dotted line) envelopes with core masses 0.9, 1.0, 1.1, 
1.2 and 1.3 M©, all of them with a similar effective temperature, ~ 4.6 x 10 5 K (~ 40 eV). Photosphcric radii are, for 
ONe50 models, between 1.9xl0 9 cm for 0.9 M© and 2.6xl0 9 cm for 1.3 M©; and for CO50 models, between 2.1xl0 9 cm 
for 0.9 M© and 2.5xl0 9 cm for 1.3 M©. Solid (ONe50) and dash-dotted (CO50) lines indicate convective zones. 
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Fig. 4. Interior structure of four envelope models with the same core mass (1.1 M Q ) and similar luminosity (4xl0 4 L Q ), 
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Fig. 10. Luminosity versus effective temperature for the high luminosity branch of CO50 and ONe50 envelopes. Time 
in days needed for the envelope to evolve between two adjacent ticks is indicated. 
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Fig. 11. Same as figure ITUI for ONe25 and ONe75 envelope models. 
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